hep-ph/0005xxx 



Neutrino Oscillations in Supersymmetry without Lepton number 

conservation and R-parity 



O ■ 

o : 

O . Tai-Fu Feng 1 ' 2 ' 3 Xue-Qian Li 1 - 2 ' 3 

(N ■ 

g ; ^CAST (World Laboratory), P.O.Box 8730, Beijing 100080, China 

Q . 2 Department of Physics, Nankai University, Tianjin 300070, China 

institute of Theoretical Physics, Academia Sinica, P.O. Box 2735, Beijing 100080, China 

(N : 
> : 

Abstract 

, With the on-shell renormalization scheme, we discuss neutrino masses up 

Cn| I to one-loop approximation in the Supersymmetry without lepton num- 

ber conservation and R-parity. Ii is shown that in this model with cx- 
. perimentally allowed parameters, Arn^ 2 , Am^ 2 and the mixing angles 

| | sin 923 1, | sin #12! which are consistent with the present observation val- 

O |' ues can be produced. We find that small neutrino mass (< 1 eV) sets a 

^ ■ loose constraint on the R-parity violation parameters in the soft breaking 

D \ terms. 
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1 Introduction 

Nonzero neutrino massQ implies that new physics beyond the Minimal Standard Model (MSM) must exist. 
During the last two decades, the minimal supersymmetric extension of the minimal standard model (MSSM)[|| 
has been studied carefully. In the supersymmetric model, the renormalizable Lagrangian can include the terms 
without lepton (L) number and/or baryon (B) number conservation. Usually, we can remove such terms by 
imposing the R-parity symmetry, with R = (—i*)3B+2L+2S j g ^ e S pj n Q £ particle). By this definition, 
the R- value of the SM particles is +1 and that of the supersymmetric partners is —1. Alternatively, we can 
dismiss the R-parity conservation and retain the B or L violation terms in the Lagrangian. Because the 
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proton decay data have set a very strict constraint on the B-violation, we presume the B-conservation and 
only keep an L-number violating interaction in the Lagrangian. In the framework of the super symmetric 
model without R-parity, some interesting phenomenological problems such as [i — ► e~f, neutrino mass and 
etc. were discussed in ||. Recently, the supersymmetric model with L-number violation has been employed 
by many authors to explain the present atmospheric neutrino experiments H |(|. The attractive point of the 
theory is that one generation neutrino can acquire mass through its mixing with gauginos and higgsinos [f?j 
and the other two- generation neutrinos can acquire their masses through loop corrections This idea urges 
people to investigate the mechanism which may induce neutrino masses ||. If in the parameter space one 
imposes a relation 62 cos 9 V + fi sin 6 V ~ where cos 6 V = v^/vl and Vd, vl are the vacuum expectation values 
of the d-type Higgs and scalar lepton, e2,/i are the R-parity conserving and violating parameters respectively, 
he can expect small v T mass due to the generalized SUSY see-saw mechanism. Recently, the authors of flOfl 
have computed the one-loop correction to neutrino masses in the MS'-scheme, they have given a very strict 
constraint on the parameter space. In this work, we carry out a complete calculation up to one-loop order 
in the on-shell renormalization scheme and investigate the constraint on the parameter space by the neutrino 
oscillation data. 

It is well known that renormalization is carried out to remove the ultraviolet divergence which appears 
in the loop calculations. At present, two renormalization schemes are often adopted in the literature. The 
first is the minimal subtraction (MS) and modified minimal subtraction (MS) schemes, these schemes are 
often applied in the QCD-calculations. The reason is that quarks and gluons are confined inside hadrons, no 
free quark or gluon can be directly observed in experiments, thus their physical masses are not measurable 
quantities. The second is the on- mass-shell subtraction scheme [jlj]], that is often used in the EW-process 
calculations. The advantage of the on-shell scheme is that all parameters have clear physical meaning and 
can be measured directly in experiments. But in general, practical calculations in this scheme are more 
complicated than in the MS'-scheme. In this paper, we would adopt the on-shell scheme for its advantage. 
Furthermore, we would perform our calculation in a strict form because any inappropriate approximation 
may impede us to understand why neutrino mass is so small, unless we have a very strong reason to take any 
approximation. Our main results can be summarized as follows: 

• Using the on-shell renormalization scheme, we find that the loop corrections to the mass matrix elements 
of neutrino-neutralino are decreasing when masses of the scalar particles turn larger. 
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• The neutrino oscillation data impose a loose constraint on the parameters of the soft breaking terms. 

An extensive analysis of the neutrino oscillations in terms of the supersymmetric model with bilinear R- 
parity violation is made in [15|. Although a different method is adopted in this work, our results qualitatively 



coincide with those of 15]. 

The paper is organized as follows: in Sec.||, we present the basic ingredients of the supersymmetric theory 
without lepton number conservation and R-parity. Using the on-shell renormalization scheme, the one-loop 
corrections to the neutrino masses are also included in the section. In Sec.|3|, we first discuss the one-loop 
corrections to the neutrino mass with only r-number violation. Then, we generalize the discussion with three- 
generation lepton number violation. We close this paper with conclusions and discussions in the last section. 
Most of the technical details are omitted in the text and then collected in the appendices. 

2 The Lagrangian and the on-shell renormalization scheme 

In the supersymmetric extension of the standard model without the lepton-number conservation, the down- 
type Higgs superfield and lepton superfield have the same gauge quantum number, we can combine them into 
a vector Lj=(Hd, Li) with J=0, 1, 2, 3. Using this notation, the superpotential can be written as 

W = fi J HLj + X JKl LjL K Ef + \ J d pq LjQ p D c q + h™H u Q p U*, (1) 

where X JKl = -X KJl . The parameters fjP, X okl and XF q are the R-parity conserving coupling constants. The 
other parameters such as //*, X kml and X k f q represent the R-parity violation coupling constants. Here we use 
the subscripts i, j, k, I, ■ ■ ■ to represent the generation indices of the leptons and quarks. 
In order to break super symmetry, the following soft terms should be introduced: 

V S oft = ^HlHu + X -L J \m\)j K L K + B J H u Lj 

+A U H U QU C + A J d LjQD c + A J e K LjL K E c + h.c. (2) 

Note, here we do not invoke the M-SUGRA scenario and absorb the superpotential parameters into the soft 
breaking parameters. In terms of W and V so ft given in Eq.(|l|) and Eq.(@), we can compute the loop corrections 
to the neutrino masses. In this paper, we will perform our analysis in the mass basis, which is independent 
of the U(4)-rotation among the superfields Lj. The relevant Feynman rules can be found in Jl3| , [l5| ] and we 
only cite those contents of the references which are necessary to the calculations of this work. 
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Now, we compute the one-loop corrections to the neutrino masses in the supersymmetric extension without 



R-parity. The general form of the self-energy for — can be written as 



(3) 



When one external leg of the self-energy is neutrino, k 2 <C rrig with mo being the mass of the heaviest 
internal-particle and we can write qj , dij , and fij as expressions of k 2 [ 12 ] : 





- r° 4- k 2 c X 




djj 


= c^. + k 2 d\j, 




e ij 






fij 


fij ^ fij' 


(4) 



Sjj's are renormalized by adding counter-terms and the renormalized T,^j EN are written as: 



T,^j EN (k) = Ey(fc) + icijirijUJ- + d*jmiOJ + + e*jjtu- + 



w 4 



(5) 



where the quantities with * are the counter parts. In the on-shell renormalization scheme they are determined 
by the mass-shell conditions 



0. 



(6) 



which yield solutions: 



d*- 



ft 



%(A0Eif*(A0lj^ = O> 

-c°- + m?4j + + i fir 

-e°- - m|c^ - m-4j - (™i + "»J)ey - rmmjflj, 

-fij ~ mimjcjj - mimjd\j - mirrijejj - (mf + m|)/^. 



(7) 



From Eq.|| and Eq.[7|, the renormalized self-energy is recast into the form 



Z% EN (k) = (mf 4 + m 2 e}j + m^/l + c^ 2 



/ O 1 01 1 1 Q \ 

+ (771^ + rnjeij + m^m^f^ + d^-fe )mju; + 

+ ( — m|c}j — mfd}j — (m 2 + m 2 )e\j — mimjflj + e}jk 2 ^k"uj- 

- miinjcjj — mimjdjj — mii7ije}j — (m 2 + m 2 )e}j + f}jk 2 ^jjtuj + 



rrij)'Eij(k)(ft - mi). 



(8) 



In the last step, we have written E^ EAr (fc) in such a way that their on-shell behavior becomes more obvious 
as 



T<ij(k) = cjjrrijUJ + + d^miU- + e^-(mjU>_ + rrijUJ + + 



(9) 



For convenience, we introduce some new symbols: 



5Z± = -mjcy - mfdjj - (mf + m|)e^ - m^f^ + ejjk 2 , 
SZij = -mimjcjj - mimjdlj - m^e^ - (m 2 + m 2 )flj + fijk 2 , 
dmfj = (mfd-j + ™l e ij + mimjflj + c^k^mj, 



(5m,- 



2 1 21 1 1 2 \ 



(10) 



Up to one-loop order, the two-point Green's function is 

Y ij {k) = ^-m t r)8 ij + ^ EN {k) 

= (f- ™f ee )<% + SZ^jku- + - tfmgw_ - <foragw H 



tree 



i.i 



+(5ij + <5Z^)(^ - m* ree - 5mg + 5Z§mf ee ')u;^ 



(11) 



where + is the renormalization multiplier for the left-handed wave function and 5ij + SZfj is the 



renormalization multiplier for the right-handed wave function. rrvf ee is the mass of the i-th generation of 



fermions at tree level. From Eq.|llj and the mass-shell conditions, we derive the loop corrections to the mass 
matrix elements as: 



5 m, 



loop 
ij 



k 2 =0 



miSZ^\ k2=m 2 + mjSZ, 



R\ 

ij I k 2 =m 2 
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= 3m* ree (mf ee ) 2 c\ + {mf ee mf ee + (mf ee ) 2 + (mf ee ) 2 )mf ™d\ 

i i'/'™i''ee\2_4ree , octree / m tree\2\ 1 , /o/ m ,tr , ee\2 OT ,tree , _tree (^tree\2\ rl fm\ 

+{{m i ) rrij +3m i [ra j ) )e ij + (3(m i ) + mi (ra,- ) (12) 

in which SZ^' R , 5mfj' R are defined in Eq.pi]. Note that the above formulae are correct for any type of fermions. 



Eq.12 is the key formulation to compute the one-loop corrections for the mass matrix of neutrino- neutralino. 



In order to compute the form factors eft , d^ 1 , e^j 1 and f®j , the one-loop self energy diagrams should be 
precisely calculated. The exchanged bosons in the diagrams can be either vector or scalar and they correspond 
to different integrals. The integral for exchanging vector-boson is 



Am Pv (k) = [ ^{iA^ JM , + m 2 + "'I (>4rW 2 )- 



dX l (2irV> (O 2 + t,(\ - x)k 2 - xrrA — (\ — x)m 2 ) 2 ^ ~ ^M^^U ~ x)Jtu ff 



(2ir) D (Q 2 + x(l - x)k 2 - xmj - (1 - x)m 2 / f 

+Dm f A^B^u (T } 

f dPQ 1 r 2x{l-x)k 2 } 

X J {2ir) D (Q 2 + xm 2 + {I - x)m 2 v ) 2 \ + Q 2 + xm) + (1 - x)m? v J 

{(2 - D)4 y )4 y )(l - x)Jku; + Dm/A^fW} (13) 



where D = 4 — 2e and /i^ represents the renormahzation scale. , B& with <r = ± are the interaction 
vertices, my represents the mass of the vector boson that appears in the loop and my is for the fermion in 
the loop. From Eq.||, Eq.|| and Eqjl^, we get 

SAmv,mt) = -iD—A^B ( L^F 2a ('mf,rnv), 
J rrij 

4(my,m/) = -iD r ^A {y) B { pF 2a {m f ,m v ), 

J m% 

e%{m v ,m f ) = -t(2 - D)A { P B^ F 2b (m f , my), 
/°(my,m/) = -*(2 - ^^B^iSjbCm/.mv), 



mi 



e|(mv,m/) = i2.A_ S_ i<3&(m/,my), 

fljimviHif) =i2A { pB ( pF 3b (m f ,m v ). (14) 



F 2a , F 2 b, F 3a and -F36 are integrals over the internal momentum of the loop and their explicit forms are given 
in appendix A. 

For exchanging scalar-boson, the amplitude is derived in a similar way and it is 



Amv(k)-(u ) 2e [^-UA^lo ) l ^±l±^ll(iB^io ) - 



= .l\ f d D Q 1 r 2x(l-x)k 2 1 

% Jo X J (2ir) D (Q 2 + xm2 + (l-x)m|)2l + Q2 + xm j + (1 _ x)m 2 / 

{4 5) ^ 5) (1 - x)JW + m f A^B^u a }, (15) 

(S) (S) 

where A a , B„ with a = ± are the interaction vertices, ms represents the mass of the scalar boson that 
appears in the loop and m/ is for the fermion in the loop. From Eq|| Eq.|] and Eq.|l5|, we obtain 

A(m s ,m f ) = i^-A ( L ) B ( L ) F 2 a(m f ,ms), 
J nij 

d°,(mg,m/) = i^-A^ ) B^ ) F 2 a(mf,ms), 
J mi 

e°(m,s,m/) = iAy B^' F 2 f, (m/ , ms ) , 

/°(m v ,m/) = ML S H S) *2&(m/,ms), 

4(™S> m /) = i— A -^ B -^ F 3a{inf,ms), 
m j 

dij(m s ,m f ) = i—*-A\_ } B y + > F 3a (m f ,m s ), 



m-i 



elj(ms,mf) = iA^ B^F 3b (mf,ms), 

fij(m v ,m f ) = iA^B^F^m^ms). (16) 
In the supersymmetric extension of the SM, the mixing of ~ originates from the following loop-diagrams. 

• The internal particles are Z ~gauge boson and neutralinos (or neutrinos through mixing) k° (a = 
l,2,---,7). 

• The internal particles are W ~gauge boson and charginos (or charged leptons) (a = 1, 2, • • • , 5). 

• The internal particles are CP-even Higgs bosons H® ((3 = 1,2, •••,5) and neutralinos (or neutrinos 
through mixing) k° (a = 1, 2, • • • , 7). 
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The internal particles are CP-odd Higgs bosons (/? = 1,2, •••,5) and neutralinos (or neutrinos 
through mixing) K° a (a = 1, 2, • • • , 7). 



The internal particles are charged Higgs bosons Ht {f5 = 1, 2, 
k~ (a = 1,2, • • • ,5). 



and charginos (or charged leptons) 



• The internal particles are up-type scalar quarks U k (a = 1,2; k=l, 2, 3.) and quark u k . 

• The internal particles are down-type scalar quarks D k (a = 1, 2; k=l, 2, 3.) and quark d k . 

The expressions of the contributions from those loop diagrams to the self-energy are presented in Appendix 
B. 



Summing over Eq.29, Eq.31, Eq.34, Eq.^T], Eq.40, Eq.42 and Eq.44 in Appendix B, we obtain the one-loop 
corrections to the neutrino-neutralino mass matrix 

Sml^oP = 6m ^°) + §m ^-) + Sm f^ + Smf>^ + 5m% + ^ + 5r^^ + 5mf > d \ (17) 

From the above analysis, we find that the one-loop corrections to the neutrino-neutralino mass matrix are 
decreasing when ms (S represents the scalar SUSY or Higgs particles) turns heavier. As a pioneer work, the 
authors of § discussed the neutrino problems in the supersymmetric model without R-parity. Even though 
our method is different from theirs, our results are qualitatively consistent with theirs, namely, as the mass 
of the SUSY particles increases, the mass of neutrinos decreases. 

3 The masses of neutrinos 

In this section, we discuss the diagonalization of the mass matrix of neutrino-neutralino with one-loop correc- 
tion. For illustrating the physics picture, we first consider a simplified model where only r-number is violated, 
and then generalize the result to the case of three-generation lepton number violation. 

3.1 A simplified model: only the r-number is violated 

In the mass basis, the one- loop correction to the neutrino-neutralino can be written as 
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(18) 



In fact, for neutralinos are much heavier than neutrinos and the loop-corrections to them are very small, so 



that here, we have ignored the loop corrections to the masses of neutralinos. From Eq.18, we can get 

1— loop 



m„ 



tree 



+ 5m. 



<5mf 5 Sm 2 ^ 5mifi 5m 2 



55 



m K o 



'25 



K 2 



'35 



K 3 



15 



m K o 



(19) 



All the off-diagonal mass terms and extra diagonal masses Srriij are given in Eq, (|l7|) . Different from the 



results in the MS'-scheme and the approach with the mass-insertion method, the on-shell scheme may give 
rise to larger corrections to the off-diagonal matrix elements 8mij(i ^ j) than to the diagonal elements for 
small neutrino masses. We will give more discussions on the difference in the last section. The corrections 
from W,Z-gauge bosons and Goldstone particles are suppressed by the small mixing between the neutrino- 
neutralino. Now, we discuss the scalar particle contributions. For convenience, we restrain our discussion in 
the basis where vl = (the vacuum expectation value of the scalar lepton). When the CP-even Higgs mass 
m H o 3> m K o , the off-diagonal correction from is 



5m i5 
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m n 



(4tt) 2 2m 2 H0 



2 ( c tu^7V" 



7 f3S 7 i(2+S) 



s wZn )(J}2 Z E Z N 
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iL E N 



c w^ E 2 9 1 



ryPl 1 . ,7/32 1 / . ry.pi 1 , 1 G 

-c w Z£ -gv 2 + s w Z^ -g vi + s w Z' E -g v 2 



4«2 r 2 



32n m 2 H0 



(20) 



As ~ 10 9 (GeV) is required, we have m H o > 3m.j. The corrections from the physical CP-odd Higgs 
bosons, Charged Higgs bosons can be discussed in a similar way, and they are suppressed by the large masses 
of the scalar bosons. For small neutrino masses, one could have another solution where all the R-parity 
violation parameters (in soft breaking terms and the superpotential) are all small, thus their contributions 
can be neglected, this scenario has been discussed by Davidson et.al fHj. 



3.2 The model with three-generation lepton number violation 

In this case, the mass matrix of neutrino-neutralino can be written as 
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(21) 



The eigenequation of the matrix is 



.3 f^tree + _ V- Sm i5 + Sm i6 + 6m i7 



1=1 n l > 
4 I W^tree 



(<5mf 6 + 5mf 7 )(m t J T ee + <5m 55 ) - 25m i5 <5mj 6 <5m 56 - 2dmi 5 5mij5m 57 



m tree 



i=l 



m o 



It is noted that even though loop is a 7th-order matrix, because of the form of the up-left corner submatrix, 
this equation can be reduced into a 3rd-order one. 

In a general theory, one can notice from the analysis given above that the lepton-number violation terms 
(bilinear or trilinear) and the corresponding soft-violating terms affect the mass matrix of neutrinos through 
mixings of neutrino-neutralino, charged lepton-chargino and Higgs-slepton. However, in this general scenario, 
there are too many parameters to reduce the prediction power of the model. Thus in our later analysis, we 
ignore the trilinear terms \ tkl LiLkEf (i, k, I = 1, 2, 3), X J f q LjQ p D^ (j,p, q = 1, 2, 3) in the superpotential and 
the corresponding soft-violating terms. 



In terms of the mass insertion method, the authors of 16] discussed the modification to the neutrino 
mass caused by influence of the R-parity violating trilinear terms of the superpotential in every details. For 
example, if only the third generation of the down-type scalar quarks is relatively light, the modification of the 
qq— loops induced by the R-parity violating trilinear terms and the corresponding soft-violating terms of the 
superpotential to the neutrino mass matrix is 

* m m 3Af Af 3 mg(^ + ^°tan/3) 

<W"(8-i)(8-i) ^2 ^2 • \ l6 > 

b 

Similarly, if the third generation of charged sleptons is the lightest, the contribution of 11— loops induced by 

J i E k 



X l ^ k LiLjE^. and the corresponding soft-violating terms is 



„- 3A* 33 A^ 3 m^f + < u°tan/?) 
*"*(8-0(8-i) §^2 ^| • ( 24 ) 

Because of the method and renormalization scheme adopted in our calculations, our results are a bit different 
from that obtained in terms of the mass-insertion method. As a matter of fact, in our calculations, if the 
R-parity violating trilinear terms are also taken into account, their effects change the mass matrix of neutrino 
through their modifications to the effective interaction vertices and thus, as a reasonable approximation, can 
be included in the effective couplings of trilinear interaction vertices. 
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Now let us turn to the numerical computations. We diagonalize the mixing matrix and obtain the neutrino 
mass eigenvalues of three generations as well as the modified neutralino masses. This is a generalized see-saw 
mechanism. 

Numerically, we adopt the basis with Vy i =0. (i=l, 2, 3) and input the soft-breaking terms as 

ml (o, 0) = -6.0 x 10 6 GeV 2 , m|(l, 1) = 6.0 x 10 6 GeV 2 , 

m|(2, 2) = 6.0 x 10 6 GeV 2 , m|(3, 3) = 6.0 x 10 6 GeV 2 , 
m 2 L (0, 1) = 0. GeV 2 , m|(0, 2) = 0. GeV 2 , m|(0, 3) = 2.0 x 10 5 GeV 2 , 
m|(l,2) = 0. GeV 2 , m|(l, 3) = 0. GeV 2 , m|(2, 3) = 0. GeV 2 . 
We also set the R-parity conserving term 

B° = 10 6 GeV 2 , 

and the R-parity violation parameters B l for three generations in the super-potential as 

B l = B 2 = 10 GeV 2 , 

and B 3 is treated as a variable which can be fixed by comparing the calculated results with data. The choice 
of the soft breaking parameters in our scenario makes the lepton-number violation effects of the 1st and 2nd 
generation fermions are much less significant than that of the 3rd generation [|17]] . 
Several other concerned parameters are taken as 

tan/3 = 5, / = 100 GeV, // = (i = l,2), /x 3 = 10~ 6 GeV, 

This choice leads to a nonzero mass less than 1 eV for the 3rd generation neutrino at the tree order. 
Then we have obtained 

m tree _ 44Q Q e y 
Xi 

m tr ee ~ 99.0 GeV, 
mK ee ~ 117.0 GeV, 

tree _ Ge y_ 
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The corresponding dependence of Am^ = \m1 T — m 2 \ and Am 2 2 = l 77 ^ ~ rn u e \ on the R-parity violation 
parameter of the third generation B 3 is shown in Fig.l. Simultaneously, we can have the mixing angles as 

| sin 9 23 1 ^ 0.26, | sin 12 \ ~ 0.034. 

In the numerical analyses, we find that the mixing angles 623, 612 do not vary much as the soft breaking 
parameter B3 changes, and this small declination can be neglected in the physical discussions. 

Therefore we can draw a rough conclusion that as the R-parity violation parameter of the third generation 
B 3 is larger than 10 5 GeV 2 , the obtained neutrino mass difference Am^ and mixing angles can meet the 
data. From Fig.l, we obtain that when B 3 > 10 5 GeV, the neutrino mass differences 

Am| 3 ~ 10" 1 eV 2 and Am 2 l2 ~ 10" 6 eV 2 . 

The recent data on the solar neutrino and atmospheric neutrino experiments have been given in The 
more precise values of the oscillation parameters at 90% c.l. are: 

Am^ T ^ (2 - 8) x 10" 3 eV 2 ; (25) 

sin 2 26^ T ~ 0.8 - 1 

and the best fit to data for the solar neutrino experiments gives several possible solutions corresponding to 
various models as 

VO : Am 2 ~ 6.5 x 10 -11 eF 2 ; sin 2 28 ~ 0.75 - 1 (26) 

SMA - MSW : Am 2 ~ 5 x 10" 6 ey 2 ; sin 2 26 ~ 5 x 10" 3 

LMA - MSW : Am 2 ~ 1.2 x 10" 5 - 3.1 x 10~ 4 eV 2 ; sin 2 26 ~ 0.58 - 1.00, 

where "VO" means "vacuum oscillation" solution, "SMA-MSW" denotes the solution of Mikheev-Smirnov- 
Wolfenstein with small mixing angle whereas " LMA-MS W" is similar to " SMA-MSW" but with larger mixing 
angle. 

Our results for Am| 3 is about 10~ 2 eV 2 for B 3 > 3 x 10 5 GeV 2 and Am 2 2 is 10~ 5 eV 2 at the scale. 
The mixing angle is sin 2 2#23 ~ 0.27 which is a bit smaller than the fitted data, but has the same order of 
magnitude. Due to the relatively larger errors in both experiments and theory, this difference is tolerable. 
The results of Am 2 2 , sin 2 2# 2 i ~ 4.6 x 10~ 3 in our calculation are consistent with the case of SMA-MSW very 
well. 
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4 Conclusion and discussion 



In this work, we employ the on-shell renormalization scheme. We find that the supersymmetric model with 
L-number violation terms can result in the expected values which are gained by fitting the solar and at- 
mospheric neutrino experiments. Even though the results are qualitatively consistent with that in the MS 
scheme, quantitatively they are a bit different from the later. To our understanding it is because the on-shell 
renormalization scheme is an over-subtraction scheme, after eliminating the divergence of the loop integra- 
tions, at any concerned order of perturbation (the one-loop order in our case), the equation of motion of 
the particle still holds, i.e. the mass-shell condition is respected. Instead, in the MS scheme, the mass-shell 
condition is violated at the concerned order of perturbation, and this may be the reason of the difference in 
the two schemes. 

When we take into account the contribution from the R-parity violating trilinear terms in our numerical 
analysis, the available parameter space would be enlarged, it is easy to find a suitable subspace of concerned 
parameters which can fit all present experimental data. In our future studies, we may extend our analysis to 
other phenomenology where the R-parity violating terms play important roles and then, maybe, the available 
parameter space is more constrained and existence of both the bilinear and trilinear terms would be necessary. 
By comparison, in this work, their effects merge into the effective couplings of the vertices for the bilinear 
interactions, so that do not manifest themselves directly. 

More studies on this subject need to be carried out, especially the parameter space for SUSY without 
R-parity is tremendously large, thus more precise measurements on neutrino oscillation as well as other rare 
B-decays may help to understand both neutrino physics and SUSY without R-parity. 
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A The definition of F 2a , F 2b , F 3a , F 3b 



The functions that appear in the text are defined as 



F 2a (mi,m 2 ) = (^) 2e J Q dx J tJ)D 



(2tt) d {Q 2 + xm\ + (1 - x)mlf 



{--1E + 4 2 > 



(47r) 2 I- e m 2 — m 2 

d D Q 1 - x 



F 26 (mi,m 2 ) = (^) 2e J Q dx J 



{2tt) d (Q 2 + xm 2 + (1 - x)m 2 ) 2 



— {- 

4vr) 2 I e 



7E + 



+ 



(m 2 — m 2 ) 2 



m 2 — 3m 2 
!(mf — m 2 ,) 

m ?( m ?-2m§)ln^+miln^]}, 



2(47r) 2 le 2(m 2 - m 2 ) 



F 3a (mi,m 2 ) = {Vwf 6 J o dx J 



d u Q 2x(l - x) 



(2ir) D (Q 2 + am? + (1 - x)m 2 



2\3 
2! 

2 ™4 „_4. 



i ; _ m 2 m 2 ln !% ^1-^2 ] 

" (4vr) 2 (m 2 -m 2 )3l in m 2 + 2 J 

x / x2, , f d D Q 2x(l-x) 2 
F 3b (mi,m 2 ) = [Hw) y o dx J — 



D (Q 2 + xm\ + (1 - x)ml) 3 

' fl 6.1 6.^4 2 24 4 2 i m l 1 

2T4 l o m i + « m 2 + n m i m 2 ~ m x m 2 - m x m 2 In \. 



" (4vr) 2 (m 2 -m 2 ) 4 l3 1 ' 6 2 ' 2 - - - - - m 2 

B The 1-loop self energy diagrams of neutralino-neutrino mixing 

• The internal particles are Z ~gauge boson and neutralinos (neutrinos through mixing) 
1, 2, • • • , 7). The coupling constants can be written as 



A a n 1 N N ~ Z^^N Z 7V h 
<±h w C w ^ S=4 

A (Z,K° a ) = _ A (Z,K° a )^ 

J 



rrif = m K o , my = m z . 
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Zn is the transformation matrix of neutrino-neutralino from the interaction basis to the mass basis. 
Corrections to the mass matrix are: 



Sm\f' K0) = \ (^m^fm^A^B^ + {{m%f e f + m^fm^ 



a=l 



tree 

j 



+ (m' K r) 2 )AL Z ' K " ) i?f' K " ) )m K oF3 a (m K o,m z ) 



{{m tr J e fm tr J e + 3miT(mi r o ee ) 2 )A^ K) B (Z ' K ° a) 



+(m^(m f J) 2 + 3(m^ cc ) 2 m^ ec )A5f'' s ° ) Bf'' s ° ) l F 36 (m K o ,m z ) 



(29) 



• The internal particles are W ~gauge boson and charginos (charged leptons) re° (a = 1, 2, ■ ■ ■ , 5). The 
coupling constants can be written as 



A (W,K a ) _ 6 f 7 tc,l 7 !,2 1 7 *a,2 yi,3\ 



V2 
i 3 

(W,K~) _ e f 7«,l 7 *j ; 2 i 1 7 a,2y*j,3\ 



1 3 

£(W,K a ) _ _e_| ^*a,l ^j,2 + _<_ ^ ^*a,2+5 ^ 4+<5 j 



= m. 



^5=0 
my = m w 



(30) 



where Z± is the transformation matrix of charged lepton-chargino from the interaction basis to the mass 
basis. The W-chargino loop corrections to the mass matrix are: 

Sm\Y' K0) = E { [Um^fm^A^B^ + ((m^) 2 + {m^f 



a=l 



+m tr J e m tr J e )A {W ^ ] 'Bf"" 1 



K - F ^{m m w ) 



+ 



tree\2^tree , i^tree ( ^tree\2\ A (W,K a ) g(W,K a ) 



+ (o(m K o J «\ + m, K o (m- K o J )A + n + J t^{m K - , mjyj | 



(31) 



• The internal particles are CP-even Higgs bosons Hg {(5 = 1, 2, • • • , 5) and neutralinos (neutrinos through 
mixing) n° a (a = 1, 2, • • • , 7). The coupling constants can be written as 



A" 



2,s w c w 
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e 


^ snn 


1s w c w 


e 


^snn 


i 2s w c w 


e 


/-<*api 
snn 


c 2s w c w 



m s = m H o, m f = m K o a (32) 

with 

= - ( 2 Z« +5 - 2ZfZ^) (33) 

5=1 

and Z^; is the mixing matrix of CP-even Higgs bosons. The CP-even Higgs-neutralino loop corrections 
to the mass matrix are 

7 5 



tree\2 



Smf' K0) = - E E {(3m*7m^i^' K " ) B^ A ° ) + ((m^ e ) 2 + (rr$* 

tree ^ree^ ^ (-^/3> K a) r)(-^fl' K a)\ j-i / \ 

+KfmX)A + B + J m K - F 3a (m K o , m^o ) 
+ (((m^ e ) 2 m£f e + 3m ^( m J ™)Vf?' K ° ) B^ A ° ) 

+(3(m^) 2 m^ + m^(m^) 2 )A^ , " a) S^ , ^y S 6(m M o,m fl o)} (34) 



The internal particles are CP-odd Higgs bosons 4o (/3 = 1, 2, ■ ■ ■ , 5) and neutralinos (neutrinos through 
mixing) k° (a = 1, 2, • • • , 7). The coupling constants are 

- ~ 'o, „ U onn 



- -i^—Conv 



with 



ms = m A o, m f = m K o a (35) 



Ct 7 = (c w Z^ - s w Z%) ( £ Zg 5 Z^ - 2Zg 2 Z^ 3 ) (36) 

5=1 
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and Zo is the mixing matrix of CP-odd Higgs bosons and neutral particles. The CP-odd Higgs-neutralino 
loop corrections to the mass matrix are 



7 5 

-3 

+ 



a=10=l V 3 1 3 

m K- F 3a{m K o m A o) 



l rrt j^Aalnf^Aa) 



+ {>s{m K o ) m K o +m K o (m K q ) )A_ B^ 



• The internal particles are charged Higgs bosons Ht {(5 = 1,2, • • • ,8) and charginos (charged leptons) 
K a ( a = 1) 2, • • • , 5). The coupling constants can be written as 



A 



(Hg,K a ) e r ,*(iaj 
— W 



' Rnk 



A 



^Lnk 



B 



e 



^Lnk 



w 1 -w 



B 



(Hp ) 
+ 



^Rnk 



S W C W 



m s = m H - 



mt = m- 



(38) 



The symbols Cink-, C^nk are defined as 
z c > f 

<) 2 

^IJK 



r a/3-r _ I \^ 7 a6\ 1 /. 7 /36 yl 2 7 08 7 yl\ 7 /31 7 7( 



7(2+5) 



E 

UK 



2>cs w c w 



7 a(5+I) 7 /3(2+J) 7 -y(4+K) V a(5+J) 7 /3(2+K) 7 7(4+J) 
Z/^ Z/_ Zj m — Zj r i Zj _ )\r 



C" 



a/37 
Rnk 



Z 



1 5 

f„ 7 */32 7 * 7 2 7 */32 7 * 7 l\ 7 */31 77 3l /X 7 <*(5+<5) 7 */3(2+<5) „ 

-j=yc w Z: + Zj n -+- s w Zj + Zj n j -f- c w Z/_|_ Z/^f j -f- v ^Su; C + 



*7l 
J 7V 



+ E 

UK 



*2es w c w 



8=3 



7 a(2+I) */3(4+ J) ,7*7(2+^) 7 a(2+J) 7 */3(4+I) 7 *~f(2+K) 
L c L N - Zj c A n 



(39) 



and Zc is the mixing matrix of charged Higgs bosons and sleptons. The Charged Higgs-chargino loop 
corrections to the mass matrix are 



5m 



(H+,k-) 



7 8 



E E \ W^m^A^^B^'^ + (m^m? + (m^ 6 ) 
q=1/3=1 1 1 3 



tree , ( ^tree\2 
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m K- F ^{m-,m H , 



+ 



^(mrf)^'" 1 ^" 1 ] F 3b (m K - , m + )}. 



(40) 



The internal particles are up-type scalar quarks (a = 1,2; k=l, 2, 3.) and quark u fc . The coupling 
constants can be written as 



A (U*,u k ) = e z *aif CwZ J2 + - Sw z jl ) - emuk Z*fZ j3 
\/2s w c w u k \ w n g w J 2y/2s w sin (3m w uk N 



A 



(U k ,u k ) 



3c w u k N 2V2 s w sin I3nn h 



ry*al r7*j3 



(U k ,u k ) _ 



B 

B m,u k ) = 



-\1< y,,-l yi\ em u k lA 

^fjk^N — TT7^ : — o ^fjk^N 

6c w u 2yJ2s w sva.f3m w u 



-Z uk [c w Z N + -s w Z N ) - ^= 



em„k 



s w sin j5m % 



fjk^N 



rrif = m u k, 



(41) 



where Zjy k denotes the mixing matrix of the k-th left-handed and right-handed up-type scalar quarks. 
The up-type scalar-quark-quark loop corrections to the mass matrix are 



ftng"» 



2 3 



i j i 3 



a=l k=l 

,treetree\ A (U k ,u k ) u {U k ,u k ) 



m K o m K o )A + .fc> + m u kl<3 a {m u k,mjj k ) 

i j J ck 



+ 



{{m tr J e fm tr J e + 3m t K 7(m7)Vf /) ^ "'^ + (m'7(m'7) 2 



+3(m^ e ) 2 <o ee )^^'" ^f"'" 'j^K*,^)} 



(42) 



The internal particles are down-type scalar quarks D k (a = 1, 2; k=l, 2, 3.) and quark d k . The coupling 
constants are written as 



A (D k ,d k ) _ e al / j2 1 7 jl\ 1 x Kkk 7 *a2 7 3( 

\/ZS w C w O z 



;(4+ JO 



3c 7 , 



r 7 *a2r 7 *jl . 1 \ ~> x E-fcfc /7*al 7 *j(4+i4") 
Z £)fc ^JV + 9 A d ^D fe ^ ' 



E"=0 
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B 



(b k ,d k ) 



+ 



"V^e r 7 a2 v i\ , 1 \ ~> v Kkh ryn I 7 



3c 



JV ' 



if=0 



5 



(D*,d fc ) 



1 „.,-i\ 1 x-^ ^Kkk 7 a2 yi^+K) 



X=0 



(43) 



where Z^ fc denotes the mixing matrix of k-th left-handed and right-handed down-type scalar quarks. 
The up type scalar quark-quark loop corrections to the mass matrix are 

2 3 



5m l?''® 



E E { [^fm^A^'^B^'^ + ((m^ ee ) 2 + (m^ e ) 2 



a=l k=l 

Jreetree^ A (D k ,d k ) „(D k ,d k ) 



-m 



tree lrec| M«„,ii > 

K o m K o )A + B + 



m d kF 3a (m d k,m £ 



+ 



tree ^2 m , tree 


3 



tree ( ^tree\2\ A (D k ,d k ) B (D k ,d k ) 



+3mX(mX)) A 



F- 
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Figure 1: Dependence of Am| on the soft R-parity violation parameter B 3 with (a) solid-line: Am 2 2 and 
(b)dash-line: Ara^. The X-axis corresponds to the B 3 (in GeV 2 ) and the Y-axis corresponds to the Am| 
(in eV 2 ). The other parameters are taken as in the text. 
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